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Abstract 

We make use of a BCS type approach based on the extended Hubbard Hamil- 
tonian to study the superconductor transition and to give a microscopic inter- 
pretation of the pressure effects on Tc in high temperature superconductors. 
This novel method suggests that the apphed pressure causes an increase of 
the superconducting gap and this effect is explored in order to explain the 
variations of Tc- Our approach is therefore beyond the scope of previously 
phenomenological models which basically postulate a pressure-induced charge 
transfer and an intrinsic term linear on the pressure. We obtain a microscopic 
interpretation of this intrinsic term and a general expansion of Tc in terms of 
the pressure. To demostrate the efficiency of the method we apply it to the 

experimental data of the Hg-based superconductors. 
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I. INTRODUCTION 



After the discovey of high temperature superconductors (HTSs) it was reahzed that the 
critical temperature (T^) could be substantially enhanced by applying high pressures. This 
very interesting effect attracted a lot of attention as it is summarized by several review 
articles which have been written on the subject The motivation for all this effort is 

not only the quest for higher temperatures but also to understand which are the parameters 
which optimise Tc besides the possibility of follow the structural changes induced by pressure 
by means of x-rays and neutron difraction techniques. It is hoped that these investigations 
lead to the preparation of new materials by chemical substitution as well as some better 
insights on these highly complex systems. 

One of the effects of the pressure which is generally accepted and well documented 
in certain materials is an increase of the carrier concentration on the Cu02 planes 

transferred from the reservoir layers. Such pressure induced charge transfer (PICT) has 
been confirmed by Hall effect and thermoeletric power measurements in several compounds 
0. Therefore this effect combined with an assumption of an intrinsic variation of (linear 
on the pressure) independent of the charge transfer was largely explored to account for 
the quantitative relation between and the pressure P and gave origin to many models 
Q-^. Some of these models also invoked that several HTSs have a versus carrier density 
n (per CUO2) diagram which satisfies a phenomenological universal parabolic curve, i.e., 
Tc = — Tjin — UopY] where Hop is the optimum n. Following along these lines we can 

write n{P) = n + An{P), Tc = Tc{n, P) and derive an expansion in powers of P, namely 

Tc{n, P) = Tc{n, 0) + ai(n)P + a2in)P^ (1) 

where the coefficients are, 

ai(n) = dTc{n)/dP - 27]T^''%n - nop)dn/dP (2a) 

and 
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a^in) = -riT^'^'^idn/dPf. (2b) 

Where the first term in ai is known as the intrinsic term and it was estimated to be con- 
stant 0-0. This approach was largely used in describing the data in the vicinity of n^p 
but fails to describe more recently measurements on different compounds with a large 
variation on the initial n values from underdoped to overdoped regime P,p!0|. Futhermore 
this phenomenological method does not provide any physical insight on the origin of the 
(charge transfer independent) intrinsic term. 



II. METHOD AND APPROACH 

We propose in this work some new ideas which are general enough to be applied to any 
family of compounds and which provide some microscopic interpretation on the effects of 
pressure. We use as starting point a recently introduced approach |[ll|] based on a BCS 
type mean field analysis which uses the extended Hubbard Hamiltonian (t-U-V) on a square 
lattice (of lattice parameter a) 

H = - ti4aCja + h.c.) + f/E^a^iT " V^ninj. (3) 

(u>.o- i {ij) 

For the sake of completeness we briefly outline the method. A gap equation at zero 
temperature is derived which has the same form as in the usual BCS theory, i.e., 

Where = —At{cos{kxa) + cosikyo)) — ji , Vgf is the Fourier transform of the potential of 
Eq. (^, which is approximately given by 

Vj:fPi U — AV {cos{kx — lx)a + cos{ky — ly)a) . (5) 



As in the BCS mean field method [|rT| the gap has the same functional form of the 
potential, namely, = Ao(cos(fcj;a) ± cos{kya))/2, where the plus sign is for the s-wave 
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and the minus sign is for d-wave channel. The chemical potential /x must be calculated self- 
consistently but as it concerns the superconducting phase boundary it suffices to approximate 
it by the value of the maximum energy (concentration dependent) in a tight-binding band. 
One then derives a finite temperature analog of Eq. and in order to determine the phase 
boundary we take T = at the limit where Aq = is applied. The zero temperature and 
the finite temperature with T = Tc gap equations are solved numerically and one matches 
Tc with Aq for a given carrier concentration n. Thus one obtains a phase diagram Tc x n 
for a given Aq which was found to reproduce well the experimental phase diagrams for 
the HTSs when the position of the attractive potential V was changed from the original 



nearest neighbor position and became an adjustable parameter |TT|. The ratio of U/V is 
not relevant to determine the phase diagram boundary as long as U ^ V but, on the 
other hand, it determines the value of Aq. The exact calculations of Aq in terms of U 
and V are not easy to be performed in a many-body system and thus it becomes a second 
parameter to be determined by comparison with the experimental Tc x n phase diagrams. 



It was shown W^ that the chosen values for Aq that reproduce the phase diagrams of the 
Lanthanum and Yttrium family of compounds are also in excelent agreement with the gap 
measurements taken from tunneling experiments and the specific heat. As concerns the ratio 
of the positions of V used for these families, it was also shown a posteriori that their values 



matches their ratio of the coherence length (as discussed in Ref. |jTT|, these are strongly 
coupled systems and the bound states are confined which is not the case of weakly coupled 
systems) thus providing a possible physical explanation for this quantity and why they are 
so different for the La and Y families of compounds. 

In connection with the above discussion we are led to propose that the effects of pressure 
are two-fold: (i)- The well accepted PICT; (ii)- The relation 2Aq = 'jksT^"'^ (7 = 3.5 
for weakly BCS and 7 fa 4.3 for the method mentioned above) suggests that if T™"^^ is a 
linear function of the pressure P (as assumed for the behavior of the intrinsic term [§|-0) 
than Aq must also be a linear function of P. As concerns the Hamiltonian of Eq.(3) this 
is equivalent to say that the structural modifications due to the applied pressure induce a 



variation on the attractive potential V (which is the most influencial parameter on the value 
of Aq). In fact the real effect of the structural changes can only be estimate by electronic 
band calculations P,|3| jr^ , p!3[] , but they are not adequate to be performed in the context of 



the strong correlations of the t-U-V Hamiltonian. 

Thus the PICT (i) implies that n{P) = n + An{P) and the assumption of a pressure 
dependent gap (ii) implies Ao(P) = Aq + AAo(P) which can be simple written as Tc = 
Tc{n{P), Ao{P)). On Fig.l we plot two curves calculated with two differents values of Aq 
to study how Tc{n) change if the pressure induces a change in Aq. Therefore to estimate 
of a compound with an initial given value of n and under a given pressure P, we perform an 
expansion of Tc{n) in terms of P. With the assumption of linear variation of n and T™"^ (or 
Aq) on the pressure, we obtain only terms up to third order, that is, 

3 

mn,P) = J2azP^ (6a) 

z=o 

with 

Furthermore one can derive simple analytical expressions for each coefficient in an ap- 
proximate way, using the universal parabolic fitting and with T^°'^{P) = 2Ao(P)/7 which 
explicits the dependency on P (assumption ii). This procedure gives an intrinsic term 
which is radically different than that used before 0-0 as well as a new third order term. 



III. COMPARISON WITH THE EXPERIMENTS 



To illustrate the entire method we will apply it to the mercury system. The mercury 
family of compounds represents a real challenge to any theory for the following reasons: (a)- 
The highest transition temperatures obtained up to date have been measured on Hgl223 
at 25-30 GPa |T^-|I^] reaching values up to 164K. (b)-The various pressure data for the 
underdoped and overdoped compounds of HgBa2Cu04_s (IIgl201) could not be interpreted 
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T0| by the models described in the introduction, (c)- The largest pressure effect on Tc with 



a change of almost 50K over a span of 20 GPa has been recently measured by one of us |T8 
in the Hg2212. 

As a preeliminary step and also in order to determine the initial parameters, we need 
to study the x n phase diagram at zero pressure. Thus we perform a calculation using 
the method of de Mello [|ll|] to obtain a Hgl201 phase diagram which is in agreement with 



the experimental results [|10l- As discussed above this method of calculation involves two 
parameters; our best result is obtained with Aq = 210K and the position of the attractive 
potential V at the 6*^-neighbor. Our results are plotted on Fig.l and, just for comparison 
and for future use, the phenomenological parabolic fitting with 77 = 50 and T^"-^ = 97K as 
used by Cao et al Thus the calculation for the Hg family phase diagram is our first 



step and it is independent of the pressure studies which we will deal next. 

To study the effects of the pressure we also plotted on Fig.l our calculation for the phase 
diagram with Aq = 2A0K. We see that near Uop = 0.16 the variations of with respect 
to Aq are almost twice as those at the extremes (low T^) and since this is the origin of the 
intrinsic term, it also varies in the same way and this behavior will be discussed below. It 
is important to notice that the two partial derivatives that appear in Eq. ( |5b| ) become two 
parameter to be determined by comparing with the Tc x P measurements for two compounds 
with different charge density n. At low pressures only the linear terms ai comes into play 
since the higher order coefficients are negligible. It is desirable (but not crucial) to start 
with n = Hop to determine dAQ/dP since at riop the charge transfer term vanishes. So with 
other set of Tc x P data at another value of n 7^ UoP, we determine dn/dP. After these two 
parameters being determined we can apply the expansion Eq. (^) to any other compound 
with different value of n. 

To illustrate the above general procedure, we will apply it to the low pressure results of 
Cao et al [|1^] for the Hgl201 and those of higher pressures of Gao et al [|T^. Our purpose 
is to show that we can describe all their results with a single choise of parameters. Thus to 
obtain the value of the two partial derivatives of Eq. (^bj) pertinent to their measurements 
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we do the following; we start with the set of data taken with the compound with n^p = 0.16, 
which has slT^x P curve at low pressures that is a straight line and from which we can infer 
that the linear term ai = 1.85K/GPa, which at Uop, is equal to the intrinsic term and from 
this we determine dAo/dP. To determine the other partial derivative we study the Tc x P 
plot made with the sample with n = 0.06. We see that the low pressure slope (ai) can be 
taken as ai = 2.6K/GPa. At n = 0.06 the intrinsic term is half of that at Uop according 
the discussion on the above paragraph and therefore, the intrinsic term is 0.9K/GPa and 
then the charge transfer must be equal to 1.7K/GPa since the sum of both terms is equal 
to «!. Using now the explicity expression for ai with the parabolic fitting with T™"^ = 97K 
and = 50 we can derive that dn/dP = 1.8 x lO^^holes/GPa. We notice that this value 
that we obtain for dn/dP is very close to others estimations P,|7|,p!8|,|9[| . Thus with this 
procedure the two derivatives which comes in the calculations of the coefficients of Eq. (|6BD 
are determined and we are set now to apply Eq. (^) to any Hg compound. Far from nop 
the charge transfer term dominates over the intrinsic one and thus the linear term ai varies 
from the 2.6K/GPa chosen above at n=0.06 up to -l.OK/GPa (at n=0.26) in the overdoped 
region. The results of our calculations for all other compounds with different densities n 
are in excellent agreement with the low pressure experimental data in both the underdoped 
and overdoped regime and are plotted on Figs. 2a and 2b. In Fig.3 we show the results for 



the high pressures measurements for the three family of mercury |jT6| at riop. As one can see 
on Fig.3 at the high pressures the quadratic and the cubic terms in the pressure expansion 
become important (for P > 20 GPa) and the agreement with the data is also remarkable. It 
is very interesting that the parameters obtained above at the low pressures yielded results 
with the maxima around 30 GPa which is the value of highest Tc for a HTSs ||16[. Our 
results can also be applied to others measurements on differents families and compounds 



of Hg which, at n^p, also yielded the maxima around 30 GPa |T^. For the first time a 
simple theory is capable to describe successfully all this ensemble of low and high pressure 
data and furthermore provides some insights on the microscopic effects of the pressure. 
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IV. CONCLUSIONS 



Thus, we conclude pointing out that our novel calculations (based on a BCS type mean 
field with the extended Hubbard Hamiltonian) demonstrate to be highly successful to de- 
scribe the effects of the pressure with just two simple assumptions. The PICT which is well 
accepted and that of the pressure induced variation of the superconductor gap which was 
introduced, to our knowledge, in this work. We hope that this assumption can be confirmed 
in the future by in situ experiments like specific heat and tunneling measurements. As this 
work came to completion, we learnt about the work of Angilella et al [TU| which also uses 
the extended Hubbard model within a BCS type approach. Furthermore they estimate the 
effects of the pressure on the attractive potential V which goes along with the lines of our 
assumption over Aq. They also obtained very good results for the effects of pressure in the 
Bi2212 family, given more support to the BCS mean field calculations with the extended 
Hubbard as a model for the the physics of the charge carriers in the HTSs. 
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FIGURES 

FIG. 1. The phase diagram for the HgBa2Cu04^s- The squares are the experimental points 
of Rcf.6. Our best fitting with Aq = 210K and the results with Aq = 240K to illustrate the effect 
of changing this parameter. 

FIG. 2. a) The lines are our caculations for the underdoped region in comparison of the 
experimental points of Ref.lO. b) The same caculations for the overdoped region in comparison of 
the experimental points of Ref.lO. 

FIG. 3. The high pressures data of Ref.l6 and our calculations (the continuous lines) described 
in the text. 
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